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Abstract 

We motivate inflationary scenarios with many scalar fields, and give a com- 
plete formulation of adiabatic and entropy perturbations. We find that if the 
potential is very flat, or if the theory has a SO{N) symmetry, the calculation 
of the fluctuation spectrum can be carried out in terms of merely two variables 
without any further assumption. We do not have to assume slow roll or SO{N) 
invariance for the background fields. We give some examples to show that, even if 
the slow roll assumption holds, the spectrum of fluctuations can be quite different 
from the case when there is a single inflaton. 



^on leave from National Taiwan University 
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1 Introduction 



As phenomenological scenarios of the early Universe, inflationary models constructed 
with a single inflaton are popular for their simplicity and consistency with existing 
observations. However, for a coherent understanding of the Universe, cosmological 
models should be embedded in fundamental particle theories. As we will argue now, 
models with many infiatons may be more natural than those with a single inflaton. 

It is often natural to have a large number of fields in various models for the early 
Universe. In Kaluza-Klein theories or string theory, it is very common to have a 
large number of moduli fields. For some PP, the enormous moduli space is the most 
important feature of string theory. For example, in the brane world scenario, there is a 
scalar field associated with each extra dimension, corresponding to brane fluctuations 
in that direction. In addition, for n D-branes on top of each other, each of these 
scalar fields is promoted to a n x n matrix. There are 6-5^ = 150 independent real 
scalar fields for 5 D3-branes in 10 dimensions, for instance. As another example, for 
the tachyon model of inflation it was pointed out that it is hard to be consistent 
with observations unless there is a large number of coincident branes pP. In fact, 
merely the 10 dimensional metric already gives 6 vectors and 15 scalars in addition to 
the metric in 4 dimensions. Recently, the possibility of inflation in pure gravity due 
to compactifications on time-varying hyperbolic spaces has attracted some attention 
and is being explored |3| . The scalars coming from the higher dimensional metric play 
the role of infiatons in these models. Based on the above examples, it may be quite 
unnatural to have a model of inflation with a single inflaton field. In fact, even the 
standard model has 3 or 4 real components from the Higgs doublet; and there are 
generically a large number of scalar fields in any grand unified theory (GUT). ^ 

Another motivation for multiple infiatons comes from the observation that in order 
to have a scale invariant spectrum of density perturbations, we need the slow roll 
^ 1 to be satisfied. This implies that we need the potential energy of 



H 



condition 

the infiatons to dominate over the kinetic energy. A naive scaling of the potential and 
kinetic energy terms with the number of infiatons indicates that this condition is 
automatically satisfied in the large N limit, as we will show in SecI21 ^ 

Inflation driven by multiple scalar fields is not a new idea. One of the examples 



""^In SU{5) GUT, there is a Higgs field in the adjoint representation with 24 components, and 
another in the fundamental representation with 5 components. There are a lot more Higgs fields for 
SO{10) GUT'S. 

^However, this scaling does not apply to assisted inflation. 
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which has been considered is the so-called "assisted inflation" scenario [HI EI, where 
exponential potentials are analyzed, and it was found that inflation is more easily 
achieved when there are a large number of scalar flelds. As another example, density 
perturbations of many scalar flelds which interact with one another only through gravity 
(but can have arbitrary self-interacting potentials) were studied in jH]. General analytic 
formulae for density perturbations in the multi-fleld inflation scenario have already been 
derived 01101111!. 

However, as we will see, it is in general a lot more complicated to deal with more 
than one fleld. In the case of a single inflaton, the slow roll condition implies that the 
slope of the potential is small. But this may not be true when there are more scalar 
flelds, as will be explained in SecE] Furthermore, for scalar flelds, one needs to solve 
N coupled 2nd order differential equations for the background flelds, and additional 
coupled equations for the perturbations of each fleld. For a large number A^ this can 
be a formidable task even with a computer. 

Despite the complexity of generic multiple inflaton models, the system can be greatly 
simplifled in special cases. In it was found that the background evolution only 
depends on a single variable if the scalars are free flelds (and the inflation is driven 
by the cosmological constant). Later it was shown 12^ that, for exponential potential, 
many scalar flelds can be assembled into an adiabatic and an entropy mode by making 
a suitable rotation in fleld space, and the same result was obtained for polynomial 
potentials 'Y2iV{(f)i) in [13J, where each fleld has a potential of the same functional 
form V. In this paper, we extend the analysis to the non-trivial case when the inflatons 
observe an SO{N) global symmetry. We flnd that only a single variable is needed 
to describe the evolution of the background geometry, and, more remarkably, only 
two equations for two independent variables are needed for the purpose of determining 
adiabatic and entropy perturbations. Note that, although we assume SO{N) symmetry 
for the interaction, the background fleld conflguration does not have to be SO{N) 
invariant. 

More speciflcally, we will focus on the case when the scalar flelds are in the fun- 
damental representation of the SO{N) symmetry. We expect that the same technique 
can be applied to more general cases with different SO{N) representations, or when 
the symmetry group is different. In those cases the number of variables needed for a 
complete description of the perturbations may be larger than two. Our experience with 
the simplest case will be valuable for those generalizations. 

Examples of scalar flelds in the fundamental representation are not hard to flnd. 
The Higgs doublet in the standard model has four real components, which transform as 
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a fundamental representation under an SO {4) symmetry, for which the Higgs potential 
is invariant. The Higgs potentials in grand unified theories do not observe any SO{N) 
symmetry in general. But if we consider the weak field approximation, and keep only the 
quadratic terms in the potential, SO{N) symmetry often emerges. For the fundamental 
and adjoint representations of SU{5), for instance, the quadratic terms in the Higgs 
potential in SU{5) GUT, Tt{WH) and Tr(S2), have 50(10) and 50(24) symmetries, 
respectively. For the brane-world scenario, there are scalar fields in the fundamental 
representation of SO{N) if there are extra (flat) dimensions. 

Imposing a symmetry on the model not only simplifies the formulation and makes 
calculations accessible, it also helps to render the model consistent with observational 
constraints on the magnitude of the entropy perturbation mode. Without the symme- 
try, each matter field gains different amounts of fluctuations according to its coupling 
to each infiaton. The predicted CMB perturbations will obtain contributions not only 
from the adiabatic mode of the metric perturbations, but also from the perturbations 
of each matter field that interacts with photons. Observational data shows that the 
adiabatic mode dominates, and thus we have to avoid excessive isocurvature perturba- 
tions. However, if the coupling of each matter field to the inflatons is SO{N) invariant, 
only SO{N) invariant operators can contribute to CMB fluctuations. Hence, we argue 
that inflationary models with multiple inflatons and with a large symmetry group are 
the most promising models to study. Similar conclusion was also reached in [T^ . 

This paper is organized as follows. After reviewing the general formulation for mul- 
tiple scalar fields coupled to gravity, we give the evolution equations for the background 
in Sec|21 and formulae for adiabatic and entropy perturbations in SecOl We studied 
perturbations for the following two circumstances: (1) theories with very flat potentials 
(SecEj), (2) theories with SO{N) symmetry (SeclHI). In the first case, if we in addi- 
tion to the condition on the potential assume the slow roll condition, we obtain the 
same scale invariant spectrum of metric perturbations as the case of a single infiaton. 
This result strengthens the case for the robustness of the connection between slow-roll 
inflation and having a scale invariant spectrum of fluctuations. 

In the second case, in which the fields have a SO{N) symmetry, we find that even if 
the background configuration is arbitrary, we still only need two variables to describe 
the adiabatic and entropy fluctuations. This is true because the coupling of the metric 
to the scalars is SO{N) invariant. 

Throughout this paper we will adopt the convention SnG = 1. 
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2 Background Equations and Large N Limit 



In this section we review the evolution of homogeneous and isotropic configurations of 
scalar fields coupled to Einstein gravity. 
For scalar fields with the Lagrangian density 

^ = E \9''''d,hduh - 1^(0) (1) 

in a spatially fiat Robertson- Walker Universe with the metric 

d.s^ = dt^ - a^{t)dx^ = a^{drf - dx^), (2) 

(the variable t denotes physical time whereas t] is conformal time) the equations of 
motion in the homogeneous background are (with an over-dot denoting the derivative 
with respect to t) 

p = K + V = 3H\ (3) 
p = K-V = -2H - 3H^, (4) 
4>i + 3H<j)j + Vi = 0. (5) 

Here, H = ^ is the Hubble rate, K = J2i ^0/ ^ ^ire the kinetic and potential 
energy densities, respectively, and Vj = 

From and (jlj), the time evolution of K and V is completely known if the scale 
factor a{t) is given 

K = -H, (6) 
V = H + 3H\ (7) 

As is well-known, not all of the equations ©, ©, and ^ are independent. Indeed, by 
multiplying (jSJ by 0/ and summing over / we see that eq. (jD) , the spatial part of the 
Einstein equation, can be derived upon using (jH)). 

Conversely, H and H are determined by K and V via © and ((7j). In particular, the 
ratio H/H'^ is determined by the ratio K/V. We have smaller H/H"^ for smaller K/V. 
In the infiationary scenario, H/H^ has to be small in order to have a scale invariant 
spectrum of density perturbation, implying that the potential energy must dominate 
over the kinetic energy. Below, we will see that potential energy domination arises 
naturally in the large A^ limit, as was already realized in jT3] (see also [T^). 

For a system of A^ scalar fields coupled to gravity, the kinetic term is proportional to 
A^. A generic potential term of the form Ylii-i ^h-i^'Ph ' ' ' 4'i„ is of order A^" because 
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it involves n consecutive sums over fields. For n > 1, which holds unless the theory is 
free, the potential energy dominates over the kinetic energy. Inflation is thus a natural 
consequence of the large limit. 

The argument above is based on the assumption that the field values 0/ and the 
coupling constants A/^.../^ do not scale with A^. The underlying philosophy is that, in a 
fundamental theory, all coupling constants are expected to be of order 1 in fundamental 
units. In string theory, for instance, the couplings among the fields are determined by 
the string coupling Qs, which is set by the expectation value of the dilaton field. If we 
include the dilaton as one of the scalar fields, then the natural assumption is that all 
couplings are of order 1. Note that one aspect of our assumption says that the initial 
values for all fields should be of order 1. This assumption is related to the types of 
arguments used to motivate the initial conditions for chaotic inflation |T2j (in models 
of chaotic inflation, the stages of inflation important for comparison with observations 
- the last 60 e-foldings of inflation - take place when the field values are of the order 1 
in Planck units). 

Notice that this large N limit is different from that of 't Hooft, where the coupling 
constants are chosen to scale with in a particular way such that quantum effects have 
a well defined limit when the classical background is set to zero. For the 0^ theory, for 
example, the coupling constant A scales like 1/A^. The target problem of 't Hooft's large 
A^ limit is very different from our consideration of the dynamics of classical background 
fields. 

Having a large number of fields means that, in general, the computations which 
have to be done in cosmology are more complicated. For instance, we will see below 
that, to calculate the metric perturbations, we need to solve A^ coupled differential 
equations. We will consider two circumstances in which the system of a large number 
of field fluctuations can be solved in terms of only two fields. 

3 Adiabatic and Entropy Perturbations 

In this section we give the general formulation for the generation and evolution of 
perturbations for models with multiple inflatons. Although equivalent formulae were 
given before in [HI QDl here we present them for completeness and also take this 
opportunity to introduce the variables n and /i that will be useful later. 
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3.1 Adiabatic Perturbation 



Fluctuations of the scalar fields 6(f)j induce fluctuations of the metric. In the longitu- 
dinal gauge (see e.g. [IH] for a survey of the theory of cosmological fluctuations), the 
metric describing scalar metric degrees of freedom is of the form 



where the two functions $ and \E' representing the fluctuations are functions of space 
and time. For matter consisting of scalar fields, there are to linear order no non- 
vanishing off-diagonal space-space components of the energy-momentum tensor. Hence, 
it follows from the off-diagonal space-space components of the Einstein equations that 
$ = \E'. Vector metric perturbations decay in an expanding Universe and will thus not 
be considered here. Tensor metric perturbations (gravitational waves) do not couple to 
matter and their evolution will hence be as in single-field inflation models. 

The Einstein equations and the equations of motion for the scalar fields at lowest 
order are given by Q and (0). The independent parts of the linear fluctuation equations 
take the form (see e.g. pUj ) 



In linear perturbation theory, each Fourier mode evolves independently. Hence, it is 
easier to solve the equations in Fourier space (i.e. after having performed a Fourier 
transformation from comoving spatial coordinates x to comoving wave numbers k). 
Note that the variables $ and 0/ have a /c-dependence which is not manifest in our 
notation. 

The Sasaki- Mukhanov ^Hj variables for the scalar fields are 





-3H^ - = li^iHi + W/) + V^(0)*, 



(9) 
(10) 




Qj = 5(1)1 + 

In terms of these variables, the above equations become (see e.g. jlOj) 



(12) 




(13) 




(14) 



(15) 



J 
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Here 



ij 



Vu- 



a" 



H 



^j(f)j + Vi(f)j + Vj(f)i 



(16) 



The comoving curvature perturbation is defined by 



As shown in [20^, in the multi-field case ( becomes 



(17) 



C 



K 



The variable ( is related to the intrinsic three-curvature on the constant energy density 
surface S via 



(3) „ -ik' 



'R 



^5 



(19) 



where '^y. = C ~ k"^"^ /{^a^H). We can use either ( or equivalently $ as a measure of 
the curvature perturbation, which represents the adiabatic mode of the metric pertur- 
bations. 

We now introduce the two gauge-invariant variables 



K 



I I 



(20) 



The curvature fiuctuation can be written in terms of these two variables as 
Conversely, we can also express k and /x in terms of $: 



(21) 



// = — -$ - {<6H^ -AH )$ - IH^ - $. 

k^ H 



(22a) 
(22b) 



In the following subsection we will see that also the entropy fiuctuation can be 
written in terms of the two variables /i and k,. 
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3.2 Entropy Perturbation 

A dimensionless gauge-invariant definition of tlie total entropy perturbation is (see e.g. 

ini) 

S = h(^J^-^-P]. (23) 
\P P J 

For N scalar fields, this yields 

2{V + 3H Zi ^^iW + 2V Zi M^h - h^) 



3(2y + 3iJE/0?)Ej0'j 
In terms of k and fi, we liave JI] 

2{H + 3HH)fi +{H + QHH) 



S 



6H{H + 3HH) 



(25) 



We see tliat botli adiabatic and entropy perturbations can be expressed in terms of tlie 
two variables k and /i. 

On the other hand, if we express the entropy perturbation in terms of $, we have 



3a^H 3HH + H [ \H 
In the long wave length limit, it turns out that [221 HHl 1211 HI] 



C ^ -^S. (27) 
P + P 



Assuming the slow roll condition {H ^ H^), the above becomes 

C ^ 3ff - I) C. (28) 

In the case of a single inflaton, S vanishes and curvature perturbations C are frozen 
outside the Hubble radius. 

In the case of multiple fields, the entropy perturbation is in general not zero, and it 
can be important in determining the magnitude of adiabatic perturbation. For example, 
as shown in ^U] and [231 5 entropy fluctuations can seed exponential growth of super- 
Hubble curvature fluctuations during inflationary reheating. In our case, when the 
infiatons are slowly rolling down the potential hill, we have p > and {p + p) > 0. 
Since H < 0, Equations and ^I7\i imply that the amplitude of super-Hubble-scale 
curvature perturbations will increase due to the coupling to the isocurvature mode. 
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3.3 Generation of Quantum Fluctuations 

Since the inflationary expansion redshifts all classical perturbations pre-existing before 
inflation, it leaves behind a quantum vacuum. In the context of inflationary cosmology, 
quantum vacuum fluctuations are believed to be the origin of the fluctuations we observe 
today. In quantizing the cosmological perturbations, it is important to identify the 
variables in terms of which the action for fluctuations takes on canonical form, i.e. in 
terms of which the fluctuation flelds have the usual kinetic term. As reviewed in [TH] . 
the canonical variables vj are related to the Qi variables via vj = a{t)Qj. Hence, 
if we perform standard canonical quantization, and expand the canonical flelds into 
creation and annihilation operators, then the expansion of the flelds Qj in terms of 
these operators takes the form (at the initial time to) 

Qiito) = (aj(to) + al(to)), (29) 



'2ka{to) 

Pj{to) = -m(to)y|(a/(to)-4(^o)), (30) 
where the Pj are the momenta conjugate to Qi 

Pjit) = a\t)Qi{t) + a\t)HQi{t). (31) 



Since the evolution of Qi ()T5j) satisfles hnear differential equations, Qi and Pj at a 
later time t must depend on the initial data linearly, 

Qiit) = ^^^^^^ J2 {(liAt)aj{to) + qij{t)a\{to)) , (32) 



Pi.j{t)aj{to) - Pij{t)a]{to)j , (33) 

where the functions qij{t)/a{t) satisfy (fT^ for all J and 

PiJ = ^-j^- (34) 

Due to (j2ni); the initial conditions are 

Quito) = Pi.j{to) = Sij. (35) 

To solve the differential equation for qijk(t), we need the initial values for both qijkito) 
and qijkito). The latter are (according to (jSH)) given by 

quito) = -ikSij. (36) 
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The two point correlation function at time t for the vacuum state defined at to, 
denoted |0)o, is 

oioiQimAtmo = ^^^\Y.^^^(^^^Kjm^'\k - k'). (3?) 

(On the left hand side, we omitted the fc-dependence of Qi and Qj, which should be 

denoted Qjk and Qjk'-) 

As we have seen in the previous section, the quantities (, $ and S, which we are 

interested in, can be expressed in terms of n and /i fl2()|) . Consequently, their correlation 

functions can also be expressed in terms of those of k, fi and their derivatives. 

In general, this calculation involves solving the equations of motion for each of the 
fields Qj. However, we will see in the following sections that, in certain cases, the 

evolution equations of k, and /i form a set of coupled differential equation by themselves. 

As a result, they can be completely determined without having to know the evolution 

of individual components Qj. 



4 Slow Roll 

The "slow roll" conditions are usually defined as 

\H\<^H\ \H\^\HH\. (38) 
It then follows from ©, © that 

r > i^, r > (39) 
For the case of a single inflaton, the slow roll conditions imply that |^| -C V^, since 

d(P^ 101 V I I ^ ^ 

Similarly one can easily show that |^^| <^ V. In this case, the slow roll conditions 
also imply that the second time derivative of the scalar field is negligible. We will now 
show that this last conclusion does not hold if there are many inflaton fields (in which 
case the use of the term "slow rolling" is in fact somewhat misleading). 

In the case of a large number of inflatons, we still need the vector 0/ to be small, 
so that K <^ V. However, we can arrange 0/, viewed as an A^-vector, to be roughly 
perpendicular to 0/ such that H = —(f)i(j)i is nearly zero. This can occur even if 0/ is not 
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small. (Note that 0/ can not be directly related to Vj if 0/ is not small.) In addition, 
since 0/ may be vanishingly small at a certain moment, Vj = V /(pi can be very large 
compared to V . Therefore, it may not be a good approximation to ignore M/j (fTHI) 
even if the slow roll conditions are satisfied, and the calculation for fluctuations will in 
general be much more complicated than in the case of a single inflaton. 

In the following we consider two situations in which the theory of perturbations is 
greatly simplified, and in which we will not have to solve coupled equations for each 
Q/, but only have two coupled equations for k and /x. 



5 Smooth Potential 

Following the previous section, we see that the "slow roll" conditions are not enough 
to guarantee the scale invariance of the spectrum of fluctuations if there are more than 
one infiatons, the reason being that the M/j terms in Equation ()16|) are not in general 
negligible. In this section we consider the case when the potential is very fiat and show 
that adiabatic and entropy perturbations can be determined by the two variables k and 
/i alone. The scale invariance of the spectrum will result if we assume that the slow 
roll conditions hold. 

Suppose that the potential is very fiat, more precisely, 

Vij < 070J, (41) 

VlJK < 0/0J0i^ (42) 

in Planck units. It follows from (PT|) that 

Mi J ~ ^ ( ^^^^^ hh + Vi^j + Vj^^ . (43) 

It is straightforward to derive the equations of motion for k, and fi 

k + 9Hk+(^!^ + 18H^ + 3if - 2|^ + I) + 2^ + ~ 0, (44) 

/i + 3ifA+ (5 + 6if+f)^ - (^ + 6if-||)/€~0. (45) 

In deriving the above we used ((Tj) and 

VjVi - -{V + 3HV), (46) 

which in turn can be derived from Equation (jj)) by multiplying it with Vj and making 
use of Equation (jUj). 
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If we also assume the slow roll conditions (jHHj) . the equations are simplified to 

k" + {e - 2re)k + 2a\fi' + = 0, (47) 

fi" + {k^ - 2n^)fi = 0, (48) 



where 



a' 

k = a^K, u = an, H = —, (49) 

a 



and the primes refer to derivatives with respect to conformal time 



The slow roll conditions imply 



dt 

dri = —. 50 
a 



H' c^H^c^ ^. (51) 



Equation ()48p is the key equation needed to demonstrate that a scale-invariant 
spectrum results, as we now show. A fiuctuation mode with wave number k crosses the 
Hubble radius at a time rjk when ~ 7i, or equivalently when 

H = (52) 

For rj <^ rjk, so that k ^ Ti, we see from Equations (j47p and (j48p that the fiuctuation 
is in the oscillatory phase. The amplitude of the canonical fiuctuation variable remains 
the same. Hence, the amplitude of the fiuctuations when they cross the Hubble radius is 
the vacuum amplitude, and as in the single infiaton case this implies that the spectrum 
will be scale invariant. 

To see this explicitly, consider that when r/ r^^, so that k <^7i, the fiuctuation k 
grows like a. The approximate solution of the equations is 

a{r]k) a{r]k) 3a{r]k) H 

where kq and /to are the initial values of k and fi prescribed at a moment tjq < rj^. 
To find the dependence on k of the two point correlation function of ( (or $), we use 
and (jSZ|). Since the dominant mode of k scales as a(t)^ (from Equation (jSBI), is 
independent of time, and hence 

where we employed (j52j) . This corresponds to a scale invariant spectrum. 
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Let us remark here that a scale invariant spectrum at the Hubble radius during 
inflation does not always imply the same for the observed spectrum of density pertur- 
bations because, as we have mentioned earlier (around ipTj) ). the entropy perturbation 
may lead to a growth of the adiabatic perturbation in time. In the present case, how- 
ever, we can solve for both the adiabatic and the entropy modes, i.e. for both /t and 
fi (both are constant on super-Hubble scales), and the joint analysis of the equations 
lead to $ having an approximately scale invariant spectrum. 



6 SO{N) Symmetry 



In this section we study another situation where we can completely determine the 
evolution of perturbations by solving for only two independent variables, without having 
to know each of the variables Qj. Consider the potential term 

V = V{B), (55) 

where 

5 = 5^02^ (56) 
I 

so that the system of infiatons has SO{N) symmetry. By multiplying Equation (0) by 
0/ and summing over J, we have 

B + 3HB + 4^ ^^1f ^ +AH = 0. (57) 
dB 

Together with Equation ((Zj), B and H form a set of coupled differential equations. 
Given B(to), B{tQ), and H{tQ), we can determine B and H uniquely. We can then solve 
for K from Equation Q. Hence, these equations completely specify the dynamics of 
the background. Since the dependence of B on H is via the functional form of V, the 
result is independent of the value of A^. 

The one special case is = 1. In this case, we have 

B = 200. (58) 

Making use of Equation ©, we obtain another equation for B and H, namely 

B^ = -8BH, (59) 

which implies Equation ()57|). In this case, we only need -B(to) and -ff(to) to specify 
the dynamics and the phase space is reduced by one dimension. Equivalently, we can 
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use Equations (|7j) and ()57p to determine B and H like in the case of multiple inflatons 
except that B(to) is now fixed by i?(to) and Hito) through Equations ((Tj) and (fS^ . 

We emphasize that although the theory itself has an SO{N) symmetry, the classical 
background is not restricted to be SO{N) invariant. Remarkably, as we will see below, 
we can still determine the metric perturbations by dealing only with two variables. 

Let us now return to the equations for the perturbations for generic N. Equation 
(|T5|l implies that 

k + 9Hk+ (^!^ + 3H + 18H^ - + f ) 

+2/^+ (6/J-2|)/i = 0, (60) 
fx+i3H- 2u)fi +(^§ + 6H+§ + u^-3Hu-u + ^B^^ /i 
-A%k - {l2H% - ±B {%Y - ^{H + 3H-){H + 6M)} n = 0. (61) 

-4log(^)^^»^™^, (62) 



Here, 



dt ^dB' 



KdBJ 



where we have used -^B = V . 

If the inverse function B = B{V) exists, then B, V, ^ and can be viewed as 
given functions of time, and thus both of the above differential equations can be used to 
solve for k and fi. In turn, these results can then be used to determine (, S and $. Our 
goal is to study under which conditions the resulting spectrum of fluctuations will be 
scale- invariant, starting with an initial vacuum state on sub-Hubble scales at the initial 
time. One obtains a scale-invariant spectrum if the Fourier modes of the canonically 
normalized flelds representing the fluctuations undergo oscillations with constant am- 
phtude while on sub-Hubble scales, and are then squeezed with amplitude proportional 
to the scale factor a{t) on super-Hubble scales (see e.g. jTSl 121] for comprehensive re- 
views). For the flelds Qj these requirements imply that we want Qj to undergo damped 
oscillations on sub-Hubble scales (with the amplitude decreasing as a~^(t), and to be 
squeezed as in the case of the previous section on super-Hubble scales). 

Assuming the slow roll conditions (jSH)), we see u and u are approximately zero and 
the Equations ()60|) . can be simplifled. In terms of the variables ()49|) and conformal 
time, these equations are 

k" + {e - 2re) k + 2a^(/x' + 2nfi) = O, (63) 
f," + fe _ 27^2 + Aa'B^) fi - Aa-'%k' + f ^ + '-^^) k = 0. (64) 
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Notice that the primes mean derivatives with respect to rj. 

A priori, these coupled equations admit the possibihty of having more than the two 
phases (sub-Hubble and super-Hubble) described above. In the case of a single inflaton, 
there is a single transition (at Hubble radius crossing) for each fluctuation mode from 
the oscillatory phase to the frozen phase. In the case of many inflatons, a fluctuation 
mode may experience more than one transition between oscillatory and frozen phases 
before the end of inflation. However, for the following explicit examples, we do not 
find multiple phases. Instead, we find that some fluctuation modes never experience a 
frozen phase, like what happens for a single scalar field with a very large mass. 

6.1 Massive Fields 

To begin, we consider the simplest case of free fields with 

V{B) = ^m^B. (65) 

In this case ^ = and = 0. The time evolution of B is determined by 

Equations ((Tj) and (j37j) . Using the definition m = m/H, Equations and can 
be simplified to yield 

k" + {e - 2re)k + 2a* ifl' + 2nfi) = O, (66) 
fi" + - 2n^)fi - 2m^n^a~'^ {k' - AHk) = 0. (67) 

Eliminating fi from the above two equations, we can obtain a fourth differential equation 
for k: 

k"" - sHk'" + {2e + 4(1 + m^)n^} k" - {se + 8(-i + 2m^)n^}nk' 

+ {k* + 4k'^n'^ + 8{l-2Th^)n*}k = 0. (68) 

Once we solve the above equation for k, we can find fi by using Equation (jUBj) . 

To get a rough idea of the solutions to these equations, we consider the following 
limits. 

(1) 7"^ ^ A; (sub-Hubble), with rh fixed. 

Equation (jUH|) can be solved order by order in Ti/k. For more detailed analysis, 
we consider the case of exponential inflation which is consistent with the slow roll 
condition. In this case 

a = -l/{Hr]), (69) 

and if is a constant. 
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(a) For m <^ 1, it reduces to the case of a smooth potentiaL We find 

k ~ CiT]^'^ cos{kri + 9i) + C2 cos^kr] + 62); (70) 
il-3H^cicos{kri + 9i). (71) 

Hence, both k and fi are in the oscillatory phase. The amplitude of the 
oscillation of /i is decreasing as a~^(t) as we need it to in order to obtain a 
scale invariant spectrum. (The relation between fi and fi is given in ()49|) .) 

(b) For m ^ 1, the fiuctuations are in the oscillatory phase with a different 
amplitude evolution 

k oc i]^^^"^ cos{kri + 61) cos(m logrj + 62)] (72) 
fi oc rj^^"^ cos{kri + 61) cos{m\ogri + 62). (73) 

{2) k (super-Hubble), with m fixed. 

Under this condition, and with (jEHl), Equation (jH^ becomes 

.„„ ^ 8 ,,„ ^ 4(1 .„ _ 8(1 ^ 8(1 - 2m^) . ^ ^ 

T] rj^ T]-^ T]^ 

Similarly, we can find the solutions of k and /t: 

l + \/9-167n^ l--\/9-16m2 

K, = C2r] + C3?7 2 + C4r7 2 ^ (75) 

^ - f 2 o -1 Z' 3 — V9 — 16m2\ 7+V9-16m2 

/X = <^ ci?7 + 3c2r] + I ^ 1 C3?7 2 



+ ( 1 1 C4r/ ^ j . (76) 

(a) For m ^ 1, the evolution equations reduce to the slow roll scenario with a 
shallow potential in Sec. El n and freeze out and the wave functions are 
squeezed as in single-field infiation. Hence, the resulting spectrum is scale 
invariant. 

(b) For m ^ 1, the large mass induces an oscillation of k and fi even on these 
super-Hubble scales: 

K oc r/"*^/^ cos(2mln?7); (77) 
/i oc r]'^''^ sin(2mlnr7). (78) 
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To summarize, in an accelerating Universe, a fluctuation mode with given k starts 
in the oscillatory phase (quantum vacuum oscillations) with k ^ Ti (Case (1)). If 
m <^ 1, then when the mode crosses the Hubble radius it freezes out and its amphtude 
grows. Thus, if m <^ 1 at the time of Hubble radius crossing, the quahtative picture of 
the evolution of fluctuations is the same as in the case of a single inflaton. 

On the other hand, if m ^ 1 at Hubble radius crossing, k, keeps oscillating even on 
super-Hubble scales. Hence, the generation of classical fluctuations is suppressed. 

Except in the case of exponential inflation, generically each fluctuation mode crosses 
the Hubble radius with a different Hubble parameter H and thus a different value of 
m. For power law inflation, i.e. a{t) oc t" with n > 1, H decreases with time, and so 
m <^ H for lower frequency modes which cross the Hubble radius earlier, and m ^ H 
for higher frequency modes. Consequently, the high frequency part of the spectrum is 
suppressed. This effect has the wrong sign to explain the suppression of the spectrum 
at large wavelengths (observed in the recent CMB data |25j) compared to what single 
fleld inflation models predict. 

6.2 0* Theory 

Finally we consider the 0^ theory, which is the interesting case in the context of the 
large N limit mentioned in Section |21 

For simplicity, we ignore the mass term in the potential and write 

V{B) = ^B'. (79) 

From ((7j) we flnd B ~ ^J^H. Making use of the deflnition A = 27A/if^, it follows 
from dnni), (El that 

k" + (fc2 - 2re) k + a* {2fi' + mfi} = 0, (80) 

fi" + + (-2 + I v^)^^} /i - |^An^a~^ {k' - lUk} = 0. (81) 

As in the previous subsection, we can combine these two equations to obtain a fourth 
order differential equation for k: 

Vln^)k"" - 8{e + -Vln^)nk"' + (2/^^ + 4k\i + -^/l)n^ + — (4v^ + ~x)n^\ k" 

3 I 3 3 J 

- |8A;^ + 8fc2(-l + 4v^)7^2 + i6(-v^ + 2A)H^}H^' (82) 

+ |a;6 + 4A;^(1 + ^ V^)7^' + Sk^l + ^A)7^^ + y (-\/I + 2A)7^^| k = Q. 
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The discussion of the growth of adiabatic perturbations parallels that of the previous 
subsection. The fluctuation starts in a phase of vacuum oscillations and then, for 
a certain region of wavelengths, freezes out at Hubble radius crossing. The other 
wavelength modes remain in the oscillatory phase even after crossing the Hubble radius. 
Let us describe each phase separately. 

(1) A; > 7^, with A fixed. 

Assuming A <^ 1, this is the short wavelength phase of vacuum oscillations in 
which k and fi evolve in exactly the same as given in the previous subsection 
(Equations ()70|) and (fTTj) ). This is far from surprising since only terms with 
leading k dependence survive in this limit. If A ^ 1, the result is the same as 
(ff^ and (ffHjl with replaced by | 

(2) k<^n, with A fixed. 

To be specific, we again consider the case of exponential inflation. In this region 
(long wavelength region). Equation (jH2I) becomes 

^ B^,. ^ 4(4 + VI) .„ ^ 12(-1 + 2Vl)^, ^ 8(-l + . ^ 

Tj Tj'^ T]'^ Tj^ 

As in the previous subsection, we can find the solutions of k and fi: 



_^ 9-16\/I 1-a/ 9-16\/I 



k = ciT] + C277 + C3?7 2 + C4?7 2 ^ (84) 
fi = HU Qcirf + ?,C2r]- + L cgr/ ^ 



,3+V9-16VA\ 7-V/9-16VI , ^ ^ 
+ I I ^ ) , (85) 



Again there are two cases: 

(a) A ^ 1. In this case, the fluctuations are frozen and their amplitude grows 
(with the leading time dependence the same, but the sub-leading terms 
slightly different from what occurs in the slow roll scenario for a flat potential 
in Sec. E}. The spectrum is still scale invariant. 
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(b) A ^ 1. In this case, the fluctuations are oscillatory even on super- Hubble 
scales, with 



This result is qualitatively similar to what was obtained in the previous subsection. 

7 Discussion 

In this paper we have studied inflationary dynamics of a system of scalar flelds 
coupled to Einstein gravity as a model for the very early Universe. In the large N 
limit, the probability for slow roll inflation increases. We showed that if the scalar flelds 
respect SO{N) symmetry, the phase space for the homogeneous, isotropic background 
is three dimensional {B,B, H), only one dimension more than the case of a single scalar 
fleld. Furthermore, for arbitrary A^, we showed that only two independent variables 
(k, fi) are needed to determine both adiabatic and entropy perturbations, and not N 
as could be naively expected. In the case of a single inflaton fleld, only one variable 
is needed since there is (on super-Hubble scales) only the adiabatic fluctuation mode. 
This simpliflcation allows us to analyze the fluctuation spectrum and flnd its deviation 
from the case of a single inflaton. Note that this model is independent of in terms of 
the appropriate variables. It will be interesting to build models of this class to explain 
experimental data and compare with single-inflaton models. 

Note that, as shown in Section El initial conditions with fleld values of order one in 
Planck units may induce slow roll inflation in the large limit. It would be interesting 
to consider the possibility of inflationary models with trivial classical background, and 
where inflation is induced purely by quantum fluctuations. In this case, a new formula- 
tion for the generation of fluctuations has to be constructed (the usual theory requires 
a non- vanishing classical background matter fleld). We leave this interesting question 
for future study (for early work on this issue see ^^)- 

Upon completion of this work, we were informed of a paper which deals with 
complementary aspects of large N inflation. 



k oc T]^'''^ cos(2A"'^/^ In rj) ; 
fl oc rj'^^'^ sin(2 In rj) . 
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